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As a simple model for single-file diffusion of hard core particles we investigate the one-dimensional 
symmetric exclusion process. We consider an open semi-infinite system where one end is coupled 
to an external reservoir of constant density p* and which initially is in an non-equilibrium state 
with bulk density po- We calculate the exact time-dependent two-point density correlation function 
Ch,i(t) = {nk{t)ni{t)) — {nk{t)){ni{t)) and the mean and variance of the integrated average net 
flux of particles N{t) — N{0) that have entered (or left) the system up to time t. We find that 
the boundary region of the semi-infinite relaxing system is in a state similar to the bulk state of a 
finite stationary system driven by a boundary gradient. The symmetric exclusion model provides 
a rare example where such behavior can be proved rigorously on the level of equal-time two-point 
correlation functions. Some implications for the relaxational dynamics of entangled polymers and 
for single- file diffusion in colloidal systems are discussed. 
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PACS: 05.70.Ln, 66.fO.Cb, 83.20.Fk, 83.70.Hq 
I. INTRODUCTION 

The simple symmetric exclusion process (SEP) |^] 
is a Markov process describing the diffusive motion of 
identical hard-core particles on a lattice. Particles hop 
randomly to nearest-neighbor sites with an exponential 
waiting-time distribution with mean tq , provided the cho- 
sen site is empty. If it is occupied, the attempt to move 
is rejected. A considerable body of exact results, par- 
ticularly for the one-dimensional case (see for re- 
views), has led to a thorough understanding of many 
fundamental properties of this model. Applications to 
specific problems include interface growth in the univer- 
sality class of the Edwards- Wilkinson equation , rep- 
tation dynamics of entangled polymers [||, or single- file 
diffusion in molecular-sized channels such as biological 
membrane channels or zeolite pores and in one- 
dimensional colloidal systems 

From a theoretical and also experimental point of view 
the main quantities of interest are (a) the equilibrium and 
stationary non-equilibrium properties of the system with 
open boundaries where particles can enter and leave 
and (b) the time evolution of the local particle density 
and their correlations in a system which relaxes after it 
has been prepared in some far-from-equilibrium initial 
state. For instance, in the study of kinetic roughening of 
an initially flat interface in a two-phase Ising system one 
obtains from the SEP exact universal scaling functions 
for the roughness [pdl . In the investigation of the repta- 
tion dynamics 1T^J13(| of an initially stretched DNA chain, 
flourescence microscopy allows for a direct measurement 
of the relaxation of the tube length which in the lat- 
tice gas approach |]l5| is proportional to the number of 
particles in the one-dimensional exclusion process with 
open boundary. Not only universal power laws, but also 
non-universal amplitudes are of interest, e.g., for estimat- 
ing the sensitivity of the coarse-graining involved in the 



lattice gas description to the microscopic properties of 
the polymer dynamics in a dense solution Pp^]. 

For a more detailed understanding of the role of 
correlations in non-equilibrium states of the symmet- 
ric exclusion process [T^ , we compute here the time- 
dependent two-point density correlation function in 
an one-dimensional semi-infinite system with one open 
boundary, connected to a reservoir of constant density 
p* . In the polymer context this equivalent to an entropic 
tensile force acting at the end segments of an entangled 
polymer chain | p^ . 

The paper is organized as follows. In order to overcome 
some of the technical difficulties connected with the cou- 
pling to boundary reservoirs we formulate the SEP in 
terms of the dynamical matrix product ansatz p"^-|2C|] 
(Sec. II). In Sec. Ill we derive a functional equation for 
the two-point density correlation function which is solved 
by a Bethe wave function. From this we derive in section 
IV the scaling form in terms of error functions. The de- 
tails of these rather involved calculations are presented in 
the appendices. In Sec. V we focus on the fluctuations in 
the total number of particles Q{t) — N{t) — N{0) which 
enter and leave the system which initially is prepared 
in an uncorrelated random state with a density pQ. We 
obtain the (expected) universal asymptotic growth law 
{Q^)-{Qf ^ A./t/^ and the non-universal amplitude 
A{pq,p*) relevant for reptation dynamics. We conclude 
with some brief remarks on the nature of the evolving 
non-equilibrium state (Sec. VI). 



II. DYNAMICAL MATRIX PRODUCT ANSATZ 

We first consider the symmetric exclusion process on 
a chain of L sites. At the boundary sites k — 1,L par- 
ticles are injected (extracted) with rates ai^L (7i,l)- In 
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terms of the transition rates Wn^n' from state n to n' the 
stochastic dynamics is described by a master equation 

^^'(21; t)= J2 [wn'~.nP{]i; t) - u.„^„.P(n; t)] (1) 

for the probabihty P(n; t) of finding, at time t, a con- 
figuration n of particles on a lattice of L sites. Here 
n = {ni, 7^2, . . . , n^} where = 0,1 are the integer- 
valued particle occupation numbers at site i. In what 
follows it is convenient to set the microscopic time unit 
To = 1. In these units the single-particle diffusion coeffi- 
cient is given by -Dq = 1/2. 

Inserting the transition rates as described above one 
readily obtains the equations of motion for m-point joint 
probabilities ( n/j^ . . . n^^ ) . These form a hierarchy of 
coupled equations where the time derivatives of the m- 
point joint probabilities are coupled to m — 1-point joint 
probabilities. In principle, this allows for a recursive so- 
lution, which, however, is difficult to obtain for m > 1. 
More importantly, unlike in a periodic system, the lack 
of translational invariance makes it difficult to obtain ex- 
act results even for m = 2 which we study in this paper. 
Hence we reformulate the dynamics in terms of the dy- 
namical matrix product ansatz p"8|-p0[| . This ansatz leads 
to a decoupling of the joint probabilities and splits the 
many-body dynamical problem into a (trivial) dynamical 
single-particle part and a (non-trivial) many-body part 
which is time-independent and which can be solved using 
the Bethe ansatz. 

The dynamical matrix product ansatz is reviewed in 
detail in ||]. Here we discuss only the main features 
relevant for the SEP. One represents an occupied (va- 
cant) site by a time-dependent matrix D (E) in a string 
DDDEDEE ... of L such matrices. The configura- 
tional probabilities P(n; t) are obtained by sandwiching 
the product of these L matrices D ot E between suit- 
ably chosen vectors ( ( | and \V)) and normalizing by 
Zl = {{W\C^\V)) where C = D + E. [Notice that 
expanding the L-th power of C automatically gives the 
sum of all unnormalized configurational probabilities and 
hence yields the correct normalization factor.] The time- 
dependent matrices satisfy algebraic relations which are 
determined by requiring the matrix product state to sat- 
isfy the master equation (|^). Expectation values of local 
observables are obtained by sandwiching suitable prod- 
ucts of the matrices D, E with C. Defining formally 
Dk — C'^^^ DC~^ one obtains for the local particle den- 
sity Pk{t) — {nk{t) ) at site k 

Pk{t)={{W\C''''DC'^-''\V))/ZL (2) 
= {{W\DkC''\V))/ZL (3) 

and for the joint probability Gk,i{t) — {nk{t)ni{t)) of 
finding particles at sites fc, I 

Gk,m^{{W\DuDiC'^\V))/ZL. (4) 



From these quantities one obtains the two-point density 
correlation function 

CuAi)^{ "fc it))-{ "fc {t)){ni{t)). (5) 

Higher-order joint probabilities are obtained analogously. 
The initial probability distribution is encoded in the ma- 
trices D(0). 

We do not review here how the matrix relations and the 
corresponding relations for the vectors \V)) are 

obtained from the master equation, but refer the reader 
to P] . The matrix C as well as the vectors ( ( W \ and 
I V Y) may be chosen to be time- independent [goj . The 
dynamical problem is then solved by introducing Fourier 
transforms 

2?p(i)=E^'''^feW- (6) 

k 

They have the simple time-dependence 

Vp{t) = e-'^^VpiO) (7) 

in terms of the initial matrix 2?p(0) and the inverse relax- 
ation times ep = 1 — cos p. These diffusive single-particle 
relaxation modes reflect the random walk nature of the 
dynamics. For calculating the local expectation values 
®li (@) it is useful to separate the static {p — 0) and 
dynamical {p ^ 0) parts in the Fourier expansion of Dk 
and to write the inverse Fourier transform in the form 

Dk{t) = il-k)V„+I+ [' ^VpiO)e-^P''~''^K (8) 

J ^TT 

In the semi-infinite system which we will be consider- 
ing, the primed integral is to be understood as a contour 
integral in the variable z = e~'^ where the contour is 
chosen in a way such that the values of the joint proba- 
bilities obey the correct initial conditions at i = 0. The 
matrices Vq and X are time-independent and yield all sta- 
tionary expectation values. The time-dependent integral 
contains the relaxational part of Dk- 

The master equation not only determines the time evo- 
lution of Dp (<) , which is given by (0) , but also requires the 
Fourier components of Dk{t), as given by (||), to satisfy 
various relations among themselves and with the vectors 
(( VF I, \V)). The dynamical part where both pi and p2 
are non-zero satisfies 

T^PiT^P2 ^ S{pi,p2)Vp^Vp^ (9) 

with 

S{P1,P2) = - ■ (10) 

The relations involving static components read 

[Vp,Vo]^0 (11) 
[Vp,I] ^2VoVp (12) 
[Vo,I]^Vl (13) 
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These relations have their origin in the bulk exclusion 
interaction between particles. 

The boundary conditions determine the action of "Dp 
on the vectors {{W \ and \ V)). One finds 

= ((W^|{Do + 2(ai+7i)X-2ai} (14) 

Q = {{W\{Vp + e^^^B,{p)V^p)] (p^O) (15) 

= {{2aL + 27i - l)I?o + 2(aL + 7l)X - 2aL} \V)) (16) 

Q = {BL{p)Vp+V_p}\V)) ip^O). (17) 



with 



2a, + 27, - 1 + e-'P 
2ai + 27i - 1 + e^P 



(18) 



Notice that in the relations (^ - (|T^) the time- 
dependence drops out. The set of relations (0), (^ - 
( p^ ) provides an alternative mathematical formulation of 
the symmetric exclusion process with open boundaries. 



III. CORRELATION FUNCTIONS 

For the equilibrium choice of boundary parameters 
a = ai = q:l = p* /2 and 7 = 71 7l = (1 - p*)/2 
the boundary relations ( p^ ) - ( p^ ) and the functions 
Bi{p) simplify considerably. With these rates, modelling 
the connection to particle reservoirs of density p*, the 
(unique) invariant measure of the process is a product 
measure with density p*, i.e., there are no density corre- 
lations between different sites. Choosing as initial state 
an uncorrelated state with density po 7^ P* leads to a 
non-trivial time evolution as the system starts to fill up 
(po < P*) or deplete (respectively po > P*)- Corre- 
lations are build up in the transient regime before the 
equilibrium state is attained. In a semi-infinite system 
this will take an infinite amount of time, and the "tran- 
sient" regime is the only relevant one. For finite sys- 
tems with L sites |l^ the system is transient for times 
t < T* (x L^. This is the regime on which we focus our 
attention. Therefore we ignore the right boundary site by 
considering the thermodynamic limit L 00. The only 
remaining length scale (besides the unit lattice constant) 
is then the dimensionless diffusion length 



L = L{t) = ^/4^/7^ro . 



(19) 



This is a dynamical length scale playing the role of a 
correlation length (see below). 

Anticipating the importance of macroscopic static ini- 
tial and equilibrium properties for the non-equilibrium 
relaxation process we introduce the basic quantities char- 
acterizing both the non-equilibrium initial state with 
density po and the asymptotic equilibrium state with den- 
sity p* . These are the density gradient 

Ap = p*- po (20) 

between bulk and boundary and the compressibility 



K = lim {{N'^)- {NY)/L 

L— *oo 



(21) 



which is readily obtained from the static two-point cor- 
relation function 



(22) 



of the uncorrelated initial and final distributions respec- 
tively. Here 5k,i is the Kronecker delta-function. Hence 



K = p[l- p) 
where p = po oi p* respectively. 



(23) 



A. Density profile and current 

The evolving density profile was computed exactly in 
Ref. @. If ai=aL= p*/2, 71 = 7l = (1 - P*)/2, one 
has 



where we use the abbreviation 



Zl 



(24) 



(25) 



in order to make the formulas more compact. Using 
equations (^ and ^ with Bi{p) = Bl{p) = e-^^P, 
one can show that {{'Dp )) obeys the functional equation 
{{Dp)) — —{{V-p)), with the solution 



{{T^p)) = E ^'^oie''"° 

feo=l 



-ipka 



(26) 



where the constants a^g are determined by the initial 
conditions. Furthermore, equations (^|) and ( [l7| ) im- 
pose constraints on the set of allowed momenta, given by 
^2ip(L+i) _ the thermodynamic limit L 00, the 

momenta p form a continuous set and this condition can 
be relaxed. For a system which is initially in an uncorre- 
lated state with density poj one has a^^ = — Ap and one 
obtains in the thermodynamic limit, after substitution of 
® in (H 



Pk{t) = Po + (p* - po)gk{t) 



(27) 



where gtit) is the lattice analog of the complementary 
error function 



9k{t) = 



hit) + 2 J2 Ip(t) 

p=k+l 



(28) 



and where Ik{t) are the modified Bessel functions ( |A.l| ). 
In terms of the scaling variable 



X — k/L 



(29) 
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the long-time behavior of the density profile is given by 
the error function, as is well-known for diffusive trans- 
port. In the vicinity of the boundary, i.e., at distances 
i <C 1 small compared to the diffusion length, the density 
profile is linear. 

Associated with the spatial variation of the density 
there is a diffusive relaxational current jk — Do(pk(t) ~ 
Pk+i{t)) which is space-independent close to the bound- 
ary to lowest leading order in time. It is convenient to 
define the current 



J = j/Do 



(30) 



in units of the single-particle diff usion coefficient Dq — 
l/(2ro). From the expansion (A. 2) of the Bessel function 
one finds 



= 2V2Ap/{t:L). 



(31) 



B. Two-point correlation function: Exact expression 

The density profile pk {t) and hence the time-dependent 
equilibrium two-point correlation functions CJl ; (t) = 
( nk{t)ni{0) ) — (p*)^ can be obtained in a straightforward 
manner from the solution of a lattice diffusion equation. 
The solution of the equations of motion for equal-time 
joint probabilities, however, and hence the calculation 
of the time-dependent two-point correlation function ^ 
is much more involved. A convenient way to circum- 
vent an explicit integration of the coupled equations is to 
make use of the algebraic representation of expectation 
values within the dynamical matrix product ansatz. Sub- 
stituting equation (^) in (^ and using the commutation 
relations (11-^, one obtains, when ai — aL = P*/2, 
71 =7L = (l-P*)/2, 

( uk {t)ni{t))=p*^+p* j ^{{-D,)) e-^P^-^"' 



dpi dp2 



2tt 2tt 

X g-HPlfc+P20-(<!pi+ep2)* 



(32) 



where {{Dp)) is given by (Pq). 

Using equations (^) and ([Tsl), we can show that 
( ( I?p2 ) ) obeys the following functional equation 



({Vp, Vp, » = -Si-pi,p2) {{Vp, V^p, )) 



(33) 



together with similar equations involving all the possi- 
ble arrangements of pi with p2 or —p2 and —pi with p2 
or —p2- Furthermore, in a finite system, these equations 
and equation ( p7| ) determine the set of allowed momenta. 



These equations are the equations obeyed by the Bethe 
wave function of a quantum spin 1/2 system with bound- 
ary fields |E1[ and they have the solution 



L 

ko<lo 



afco,io*Pi,P2(^0, ^o) 



(34) 



where the constants aka,ia ^''^ determined by the initial 
conditions and ^P^j^^p^ (fco, ^o) is the Bethe wave function, 
which is given by 

^p„p,(hM) = e^P^'^^+'P-'" + S{pi,p2)e'P'''°+'P''" (35) 

_ g-ipiko+ip2lo _ S{—pi P2) e^P^^°~^P^''° 

- S{pi,p2) e~'P^^°+'P'^'> 
+ S{~pi,p2)e-'P'''°-'P''« 

- S{-pi,P2) S{pi,P2) 



^ipiko-ip^lo 



+ S{-puP2)Sipi,p2)e-'P'''°-'P''° 



Using equation (|2j), one can write the two-point corre- 
lation function as 



dpi dp2 



{{{T^P.T^P.)) 



2tt 2tt 

{{Vp, )){{Vp, )) )e-'(fi'=+P2')-(^pi+<^P2)* 



(36) 



CkAt) = 



where {{Vp^ Vp^ )) is given by (|3^) and (( )) by (|2e 

If we now choose the contour of integration such that 
at t = the integral over ^fp^^p^ (fep, ^0) in ( ^ is equal 
to 5k,ka ^i,ioi then the condition that the initial state is 
uncorrelated Ck,iiO) = 0, gives ako,io = afco "-lo = (^P)^- 
Substituting this result in ( |3^ ) and using this equation, 
together with ( p6| ) in (|3^), we obtain for the semi- infinite 
system, when L ^ 00, 



CkAt) = (Apf {e-^Ali+ki2t) - Ii^ki2t) ) 



(37) 



. 2n 2n 

ko<lo 

X [{SipuP2) ~ I) e'P'''°+'P'^" 

-iSi-pi,p2)-l)e'P'''"-'P''" 

- {S{pi,P2) - l)e-*P^'=o+*Pi'" 

+ {S{-pi,p2) - 1) e^'P^'^'-'P'^" 

-iSi-pi,p2)S{p,,p2)^l)e'P'''°-^P^^'> 

+ {S{-P1,P2) S{pi,P2) - 1) e'^P^>'''-^P-'» ] } , 

where the first term comes from the summation terms of 
((I?pj ))((I?p2 )) with fco = ^0- We remark that for non- 
interacting particles (no hard-core repulsion and hence 
no exclusion, but otherwise identical hopping dynamics) 
the correlation function has only this term, but with a 
different amplitude oq = ( )q — Pq — po determined by 
the initial distribution only. The exclusion interaction 
gives rise to the double sum over fco, lo with the terms of 
the form S'-l, SS' ~1. 
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To extract more detailed information from the exact 
expression (|3^), one needs to perform the double sum 
over ko,lo and to determine the contour of integration 
such that ( p7|) obeys the initial condition. This is shown 
in Appendix M. As a result, we obtain the more compact 
expression 



where 



Ck,i{t) = -{Apf [Fk+i-iA- 
—Fi-k-i,k+i{t) — 



k{t) + Fk+Ll-k+l{t) 
Fl-k,k+l+l{t)] (38) 



with 



2x 1-7^/2 



F™„(a;) = ^^ / ' d^cos^ / dve^^'°''^x 
I" Jo Ja 
Im{'2{x - v) cos e) COS [v sin (261) + n9] , (39) 

which will be used in the next section to obtain the scal- 
ing behaviour of the two-point correlation function and in 
section ^ to determine the particle number fluctuations. 



IV. SCALING BEHAVIOUR OF THE 
TWO-POINT CORRELATION FUNCTION 

The expression ( ^8|) is our starting point for analyzing 
the particle number fluctuations (see next section). For 
investigating local correlations it is more convenient to 
write ( |3^ ) in another way 

Ck,iit) = -(Ap)2 [Pkj-iit) + FkAt) 

-F^kj~i{t)-F-kA{t)\ (40) 

with 



(1 - «)"-™+^/,„(t(l - v))I,M^ - v)). (41) 

One can sho w th at ( ^ ) holds if we use the integral 
representation ( |A.l|) of the modified Bessel functions ap- 
pearing in the definition of Fm,n{x), as given by (39). 
We can then perform the integral over v. Representing 
the integral over ^ as a complex integral over an appro- 
priate contour we then obtain, after performing some 
expansions and using the identities (A. 14) and (A.S), 
i)/2jv™+n-i)/2(a;) which, when sub- 



[Jni-n+l) /2Jm+n- 

stituted in (^, yields gg). 

For t S> 1 the main contribution to the integral comes 
from small values of w. In order to obtain the scaling 
behavior we define in analogy to ( p9|) a second scaling 
variable y — l/L and substitute the integration vari- 
able u — vL. To leading order in time one then has 
(1 — v)''~^ = exp [—u{x — y)]. With (A. 2) we obtain the 



scaling form of the correlation function 



C^yit) = -{Apyj-[Rii,y) - R{-S:,y)] (42) 



R{x,y) 



^erfc((2/ x)/y/TT). 



(43) 



and erfc(.) is the complementary error function (Fig. 1). 
As time increases, the range of correlations increases pro- 
portionally to y/t, but the amplitude decreases in the 
same manner. The negative sign of the correlation func- 
tion signals anticorrelations typical for the exclusion ef- 
fect 0. In a finite system with a fixed density gradi- 
ent between the two boundaries imposed by the coupling 
to two different particle reservoirs these anticorrelations 
persist |l^ . They extend over the whole lattice and have 
an amplitude inversely proportional to the system size. 




FIG. 1. Scaling part R = R{x,y) — R{—x,y) of the 
two-point correlation function as a function of the scaled bulk 
coordinate x and the scaled lattice distance f — y — x (full 
lines) . For visualization purposes also the plane R = (equi- 
librium case) is shown (broken lines). 

In the vicinity of the boundary {0 < x < y <ti 1) 
one has R{x,y) — R{—x,y) — 4i(l — ?/)/7r. In terms of 
the current (^Tj) and the diffusion length L{t) ( p^ ) the 
boundary correlation function is thus given by 



C = -Lfi{l - y). 



(44) 



Corrections, which can easily be obtained from the exact 
scaling function (^), are of third order in the scaling 
variables. 

In order to disentangle the effects caused by the initial 
distribution and the exclusion interaction, a comparison 
with non-interacting particles is again instructive. The 
scaling form of the correlation function is similar in struc- 
ture, but one has R{x,y) = e~*^^~^-' and a different 



amplitude oq / %/ Ant (see remark after Eq. (37)). For the 
same uncorrelated initial distribution as considered for 
the exclusion process one has Oq = — Pg < 0. Remark- 
ably, anticorrelations develop even though there is no ex- 
clusion. However, the amplitude is different and, unlike 
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in the exclusion process, these anticorrelations vanish in 
a finite system driven by a boundary gradient as time 
tends to infinity. 



V. PARTICLE NUMBER FLUCTUATIONS 

We define as Q{t) = N{t) - N{0) the net number of 
particles that have entered or left the s yste m until ti me t . 
The mean {Q{t)) is evaluated using (^), ( |A.4[ ), ( [A^ ) 
and one gets 



[Q{t))^^pY,gk{t) 

oo 

= Ape-*^(2fc-l)4(i) 



= Ap 



fc=i 



te-*(/o(t)+/i(<))- 



(45) 



The mean grows asymptotically with the power law 



V TT 

characteristic for diffusive processes. 
The variance 

a^{t)^{QHt))-{Qit)f 

may be split into three different parts 

a'it) ^ {N^it))- {N{t)f + {NHO))- (N{0)f 
-2i{Nit)N{0))~{N{t)){N{0))). 



(46) 



(47) 



(48) 



Since the initial state is a product state one has 
{N{t)NiO)) - {Nit)){NiO)) = no^iNit)) where 
(see (p3|)) is the compressibility of the system in the ini- 
tial state. 

On the other hand, because of the exclusion principle, 
one can write 



( N^t) ) - ( N{t) f = ~{ApfKit) + J2{pkit) - plit)) 

k=l 

(49) 

with the double sum 

oo oo 

Using 



k=l l=k+l 



(50) 



— {nk{t)) = l-gk{t) 
dpo 



(51) 



we rewrite (Eq) in the form 



<^Ht) = + + i^Pf)9kit) - {Apf gl{t)] 



k=l 



-{Apfm 



(52) 



convenient for studying its asymptotic behavior. The 
sum over gk has been calculated above (^5|). The eval- 
uation of the other sums is rather technical, the details 
are presented in Appendix |^. One finds the following 
asymptotics 



K{t) - (3- 2\/2)i 



2-1) 



fc=i 

and therefore 



a''{t)= A{pa,p*)Vt 



with 



Ko + K* + (3-2V2)(Ap)' 



(53) 
(54) 

(55) 
(56) 



The amplitude A is symmetric under interchange of the 
macroscopic quantities po,p* and convex in the physical 
domain < pq, p* , < 1 with a local maximum at po = 
p* = 1/2 (Fig. H). For an initially completely filled lattice 
(po = 1) we recover the result presented previously [ p^ . 
For the trivial case of non-interacting particles one has 
(\/2 — l)po instead of the gradient term in (|5F 




FIG. 2. Amplitude A of the particle number fluctuations 
as a function of the boundary density p* and the initial bulk 
density po. 



Within the Rubinstein model for reptation |15| and 
its extension the particle number of the symmetric 
exclusion process is proportional to the tube length of a 
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polymer and hence proportional to the experimentally ac- 
cessible visual length of a flourescence marked entangled 
macromolecule such as DNA. An non-equilibrium initial 
density po > P* corresponds to a stretched conforma- 
tion which may be approximated by dragging a molecule 
through a dense solution with optical tweezers 10. It 
has been shown that the relaxing tube length calculated 
from the symmetric exclusion process is in good agree- 
ment with experimental data in the universal initial-time 
regime The particle number fluctuations (the lattice 
sum over the two-point density correlation function) that 
we obtain yield the evolution of the tube length fluctua- 
tions of the polymer chain. For an initially fully stretched 
polymer (po ~ 1) the result (|5|), ( p6| ) has been discussed 
in a recent publication . We see here that a partially 
stretched chain displays qualitatively similar relaxational 
behavior. In contrast to the Langevin approach used in 
Rouse-based standard reptation theory ||l^ where A de- 
pends solely on the equilibrium tube length fluctuations, 
the exclusion model predicts a dependence also on the ini- 
tial state (via the compressibility kq) and the amount of 
stretching Ap. These features arise from the static inter- 
actions between particles (hard-core, in our case) which 
are neglected in the purely entropic Rouse model. 



VI. CONCLUSIONS 

The main results of our study are the exact expres- 
sion (40) and scaling form ( ^2| ) respectively of the time- 
dependent two-point density correlation function and the 
asymptotic variance ( p5| ) in the number of particles that 
have left or entered the system up to time t. It turns out 
that the non-equilibrium behavior of the model is largely 
determined by the dynamical diffusion length L and by 
three static macroscopic quantities, viz. the compress- 
ibilities Ko,K* of the initial and equilibrium states resp., 
and the density gradient Ap between the (equilibrium) 
boundary and the (initial) bulk density. With regard to 
polymer re pta tion this result ( |56| ) supports our previous 
conclusion |lq ] that standard reptation theory is an over- 
simplified model of the relaxation process of stretched, 
entangled polymers. 

On a local level we find an inverse relationship be- 
tween the range ^ cx L of (anti-) correlations and their 
strength (x l/L. It is interesting to quantitatively com- 
pare the expression ( ^ ) for the correlation function in 
the boundary region of the system with the stationary 
non-equilibrium correlation function C* of a finite system 
of L sites with two different reservoir densities p~ = p* , 
p+ = po as considered in the Section II. The external 
density gradient imposes a stationary current which is 
given by J* = Ap/L. In terms of the scaling variables 
X = k/L, y^l/L one finds C* = -L{]*fx{l-y) @j20| 
which is of the same form as (^) . This result suggests 
that the state of the open system in the boundary region 
(i.e., at distances small compared to the diffusion length 



L) is similar to the non-equilibrium steady state of a fi- 
nite system of size L = L. Hence we can identify three 
distinct length scales where the system displays different 
behavior. On the scale of the lattice constant a = 1 the 
system (in the boundary region) is in local equilibrium 
as is the bulk state of the finite stationary system . 
On intermediate scales a <C r <C i the system is locally 
(i.e., in the boundary region) stationary, but not in equi- 
librium. On large scales r ^ L the system is neither 
in equilibrium nor stationary, but displays relaxational 
behavior and dynamical scaling. 

It is no surprise that the qualitative features of the 
relaxation process in simple symmetric exclusion can be 
described in terms of dynamical scaling with the diffusion 
length L{t) and the universal power law characteris- 
tic for diffusive dynamics. Yet it is gratifying to have a 
simple, but non-trivial model, where not only scaling the- 
ories can be verified explicitly, but also scaling functions 
can be calculated. An interesting open problem remains 
the question whether these results can be obtained from 
more widely applicable coarse-grained hydrodynamic ap- 
proaches for stochastic interacting particle systems. 
A direct experimental study of the questions addressed 
here appears to be feasible by studying colloidal parti- 
cles in a set-up similar to that used in ||^. In such an 
experiment the influence of direct particle interactions in 
addition to pure hard-core repulsion can be studied. 
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APPENDIX A: MODIFIED BESSEL FUNCTIONS 
AND ELLIPTIC INTEGRALS 

1. Modified Bessel functions 

Here we list some useful properties (see e.g. |2^) of the 
modified Bessel functions 



/„(t) 



1 

2^ 



i(pn-\-cos (pt 



(A.l) 



with integer index n. 

(i) Asymptotic behavior {t oo, u — /t finite): 



e~*/„(t) 



(ii) Recursion relations: 



1 



V(2i) 



(A.2) 
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In{t)=I-n{t) 
2nI^{t)=t{In-l{t)~In+l{t)) 

2 ^ Init) = In-lit) + In+l{t) 

at 

(iii) Summation formulae: 

oo 
n— — oo 

oo 

E In{t)ln+m{t) = I„,{2t) 



(A.3) 

(A.4) 
(A.5) 



(A.6) 
(A.7) 



n— — oo 



(iv) Integrals: 

For m, 71 integers, one has 

Im{t)In{t) ^ - / (i6' cos((n =F '7i)6')/„±„(2t COS ( 
7^ Jo 



One defines ||2^ the functions fq.s{t) where g is a pos- 
itive integer and s is an integer as 



US) ^ E 

p=0 



p + q-1 
P 



2''(-irW2pW- (A.9) 



One can also show, using the integral representation of 
the modified Bessel functions, that for s > 0, one has 



Usit) = 



1 



du Is{u) {t — u) 



9-1 



(9-1)! J, 
The following useful equalities also hold 

UAt) = \{f,+ l.s + l{t) + fq+l,s-l(t)) ., 
fq+l,l{t) fq.n{t) T , 



(A.IO) 

(A.ll) 
(A.12) 

(A.13) 



where the first equality follows from the integral repre- 
sentation of .fo.s(t) , the second follows from integration 
by parts of ( A.IO ) and the third follows from the two 
above. 

One can also show |2^, that the following identity 
holds 



oo „i 

E Ip+mix)yP (1 - w)™ e-*^«^'+^«^ 

X I^_^{x{l - v)) , (A.14) 
where m is a positive integer. 



2. Elliptic integrals 

The following relations are used in the calculation 
of Laplace transforms: 



r/2 



d6- 



(1 - a2 sin^ 61) \/ 1 - sin^ ( 

7r(l- Ao(<^, k)) 



2^a2(l-a2)(a2-fc2) 



(A.15) 



where k < a, sm(j> = — c?)j (1 — fc2) and Ao((/', A:) is 
given in terms of elliptic functions by 

Ao(0,fc) = -{E{k)F{(j),k') + K{k)E{(t),k') 

-K{k)F{^,k')), (A.16) 

where K{k) and E{k) are the complete elliptic integrals 
and F{(f), k'), E{4>, k') are the elliptic integrals of first and 



(A. 8) second kind, with k' = vT— fc^- 



APPENDIX B: DERIVATION OF THE EXACT 
EXPRESSION FOR THE TWO-POINT 
CORRELATION FUNCTION 

In order to derive the exact expression for the 
two-point correlation function from equation (37), one 
needs, as stated above, to perform the double sum 
X]fe^=i X]i^=fco+i ^ — ^ ^^'^ then to determine the con- 
tour of integration of the double integral in this equa- 
tion which will yield the correct initial condition, namely 
Cfc,i(0) = 0. In order to perform the first step, the 
key point is to realize that one can write the factors 
5(^1,^2) - 1, S{-puP2) - 1 and S{-pi,p2) 5(pi,p2) - 1 
which appear in ( |37| ) in a way such that the sums over 
Iq coming from each of these terms can be written as 
the difference of two sums starting at neighbouring argu- 
ments, e.g. Iq = ko + I and ^0 = ^0 + 2, and can thus be 
easily performed using the telescopic property of sums. 
Furthermore, after some tedious but straightforward al- 
gebraic manipulations one can show that the sums over 
fco can also be performed in the same way, i.e. using the 
telescopic property. When performing this second sum, 
one also generates one extra term which cancels exactly 
the first term of (|37|). One obtains, after interchanging 
Pi and p2 , the following result 

CAt) - -2(Ap)^ /' /' |i |^e-(^-+^-)* (B.l) 



1 



1 - 2e~'Pi -I- e-*Pi-'P2 



-ip^k — ip-i I 



1 



-ip2k—ipil 



1 — 2e^*Pi + e-*Pi+'*P2 
Now we need to determine the appropriate contour of 



integration in (B.l). But in fact, one does not need to de- 
termine it explicitly. Assuming that the contour includes 



8 



the origin, one can use the identity 1 j x = L"" da e to 



represent each of the denominators of (BJ) as an integral 
over a and then formally expand the resulting exponen- 
tials under the integration sign in powers of e"'^^ and 



for the first denominator, 



and 



-ipi+ip2 



for the second denominator. If one then performs the in- 
tegrals over pi, p2 and a, one obtains the following result 
for 



J5 + 9 
V 



29+1 (3 2) 



p,(j=0 

X ( /p+fe {t) Ip+q+l (t) + Ip+k [t] Ip+q+l+l (t) 
- Ip-k{t) Ip+q+l{t) - Ip-k{t) Ip+q+l+l{t) ) . 

Notice that the existence of the expansions depends on 
the convergence of the resulting series, which implicitly 
fixes the contour. It can be easily checked from the prop- 
erties of the modified Bessel functions that this expres- 
sion does indeed obey the initial condition Cfc_;(0) = 
(notice that k < I). 

This exp ression is still rather cumbersome to use. If we 
apply (A. 8) to products of two modified Bessel functions, 
we will obtain, with ^1 = 1 — k, n = l-\-k^r = q + 2p 



Ck.iit) 



2(Ap)2 _ 



TT 



/.7r/2 °° / , \ 



p,q=0 

^r+i+{n~i) ( 2i cos 6 ) cos( {q + m)9 ) 



X (/^ 

-t- /r+i+n( 2tcos9) cos( {q + 1 + m)9 ) 
- 4+i+(m-i)( 2^008 6*) cos.{{q + n)9) 
-Ir+i+mi2tcos9) cos{{q + l + n)9)) . (B.3) 



One can now write (B.3) in terms of the fq,s{t) func- 
tions which were defined above. One has 



Cfe,i(t) = - 



d9 



(B.4) 



[ fq+i,i+k-i{2t cos 61) cos( {q + l-k)9) 

q=a 

+ fq+id+k (2i cos 9) cos( {q + l-k + l)e) 

- /g+i j-fc-i (2i cos 9) cos( {q + l + k)9) 

- fq+U-k{'^t COB 6) COB{{q + l + k+l)9)] 



Using the integral representation ( |A.10 ) and summing 
over q one obtains, after the substitution u = 2{t—v) cos9 
in the integral over u, the solution ( |3^ ) where the func- 
tions F,„_„(x) are given by (^9|). 

Since the steps which led from (B.l) to (p3.2p are only 
formal, one should check explicitly that (|38D is indeed a 
solution of the equations of motion for the joint probabili- 
ties ( nfc(t)ni(i) ). This is trivial for / 7^ fc-|-l. ForZ = fc-|-l 
one can show that the unphysical amplitudes Ck,kit), 
Ck+i,k+i{t) which are generated by the time derivative 
of {nk(t)ni^^i{t) ), obey the following identity 



Ck,k{t) + Ck+i,k+i{t) — 2CkM+i{t) = 
-{Apfe-^' {Ik{t) + Ik+i{t)f 



(B.5) 



which cancels exactly the unphysical contribution com- 
ing from the term {pk{t) — pk+i{t))'^ (see equations 
( ^ ) and (|2^)) which also appears in the equation for 
(nk{t)nk+i{t)) , thus showing that {nk{t)nk+i{t)) obeys 
the correct equation of motion. This identity ca n be 
proved by considering the expression for the Ihs of ( B.5 ) 
as given in terms of the integral representation ( |3^ ) . Af- 
ter some cancelations between the terms, one uses the 
identity 



dv 



(gi,(i+cos(2e)) cos(i;sin(26l) + sO) ) 



2 21.(1+008(29)) cos(wsin(26') + (s + 1)6*) cos 6* (B.6) 

to integrate the r esulting expression by parts. Applying 
the identity ( |A.8|) to the remaining integrals, in order to 
transform them into products of modified Bessel func- 
tions yields the desired result. An altern ative route to 
derive (|B.5|) is to us e the representation (B.4) and the 



equalities ( |A.ll ) to ( A. 13 ) 



APPENDIX C: ASYMPTOTICS OF cr^ 
1. Sum over gi(t) 

This expression arises in the summation over ( nk{t) )^ 
which forms part of the 'dynamical compressibility' 
{N'^[t)) — {N{t))'^ entering the expression for cr^(i). It 
is convenient to split this sum into three different parts. 
One has 



k=l 



k=l \fc=l / 

oo oo oo 

k—l p—k q—k 



The first part can be e valuated using ( [A.3| ), (A. 7), the 
second part using (A.6). To evaluate the third part we 
rewrite the summation 



c« oo oo 



k-1 



EEE^f w^^w = E [f'^'it) + 2j2pipm{t) 

k=lp=kq=k k=l \ p=0 ) 

(C.2) 



and apply (A.3), (A.4), ([Aj|, (A.7). One obtains 



OO OO OO 



EEE^fW^«(*) = 2 [e*(^o(t) +/i(0) 

h—\ p—k q—k 

-(/o(2t) + /i(2i))]. (C.3) 



Putting everything together yields 
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Y^glit) = 2t [e-*(/o(<) + hit)) ~ e-'\lo{2t) + h{2t))] 
-(1 - 2e~'lo{t) + e-2*/o(2t))/2. (C.4) 



where Q!^(s) 



k=l 



The second part of this expression is subleading in time 
and can be ignored in the study of the asymptotic be- 
havior. 



2. Laplace transform of K{t) 

In order to perform the sum on the right hand side 
of (50), we consider the Laplace transform of the cor- 
relation function Ck,i{t) as given by (p8|). Since the 
expression ( ^8| ) involves a convolution of two functions 

/o /(^ ^ ^)9W)i such transformation simplifies consid- 
erably the calculations, because the Laplace transform of 
such a convolution is the product of the Laplace trans- 
forms of the two functions. Writing the Laplace trans- 
form as Ck,i{s), one has 



Ck,i{s) 



4(Ap)2 



Tr/2 



cos 9 



^(s + 2)2 -4cos2 6' 
1 



((s + l)2-2(s + l)cos(26l) + l) 
X [^k+i-i.i-k{s) + ^k+i,i-k+i{s) 

~ ^l-k-l,k+l{,s) — Ti-k.k+1 + li.s)] , 



(C.5) 



where the functions !Frn,n{s) are defined by 



+ 2- ^(s + 2)2 -4cos2 



I 2cos6l / 

X [(s + 1) cos(n6l) - cos((n - 2)61)] . (C.6) 



An additional advantage of ( |C.5D with respect to (|3§| ) 
is that the sums over k and I now reduce to summing two 
geometric series, due to the form of Tm,n[s). Performing 
such sums, one obtains after computing some standard 
integrals using the residue theorem, the following expres- 
sion for the Laplace transform K{s) of K{t) 



4(3s + 4) 
^ 7r(s + 2)3 

X / de 

Jo 



(C.7) 



sm 



(1 - a2(s) sin^ 61)^/1 - k^{s)sin^9 



7r(s + 2)2 
X / d9 



sm 



(1 - a'2(s)sin2 9)^1 - k^{s) sin^ 6* ' 



4(s-|-l) 
(s-|-2)2 



a'2(s) 



sT2 and fc2(s) = 



Since the last two terms are of the form ( A.15[) , one can 
now expand the elliptic integrals at small s ||2^. The 
most singular terms of this expansion, i.e. the terms 
which diverge as 1/5^/2 at small s diverge like Vi at large 
t in the time-domain, i.e. after inverting the Laplace 
transformation. Collecting all these leading order terms 
in the expansion of K{s) yields the result given in equa- 
tion dH). 
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